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We develop the formalism necessary to study four-point functions of the cosmic microwave back- 
ground (CMB) temperature and polarization fields. We determine the general form of CMB trispec- 
tra, with the constraints imposed by the assumption of statistical isotropy of the CMB fields, and 
derive expressions for their estimators, as well as their Gaussian noise properties. We apply these 
techniques to initial non-Gaussianity of a form motivated by inflationary models. Due to the large 
number of four-point configurations, the sensitivity of the trispectra to initial non-Gaussianity ap- 
proaches that of the temperature bispectrum at high multipole moment. These trispectra techniques 
will also be useful in the study of secondary anisotropies induced for example by the gravitational 
lensing of the CMB by the large scale structure of the universe. 



I. INTRODUCTION 

Beyond the power spectra of the cosmic microwave background (CMB) lies the relatively unexplored territory of 
non-Gaussian statistics. Studies of its non-Gaussianity hold the potential to reveal physics at the two ends of time. 
Non-Gaussianity in the primary anisotropies from the recombination epoch can test the inflationary model of the 
origin of fluctuations (e.g., EL S, 0, 0, Non-Gaussianity in the secondary anisotropies, arising during the transit 
of a CMB photon through the large-scale structure of the universe, probes the nature of the dark energy and dark 
matter (e.g., | J J, | 0). 

The primary challenge facing non-Gaussian studies of the CMB is the selection of an appropriate statistics. The 
term "non-Gaussianity" tells us what the distribution is not, not what it is. Like the power spectra, the higher-point 
correlations of the multipole moments of CMB fields provide a set of statistics with definitive predictions in the cases 
of cosmological interest. Unlike the power spectra, there are a large number of potential observables, associated 
with the distinct configurations of the points, requiring the development of new techniques for their prediction and 
estimation. In particular, it is important to identify the symmetry properties of the spectra to build optimal statistics 
for the detection of non-Gaussianity. 

Non-Gaussian signatures in the three-point function or bispectrum of the temperature distribution ||, |ll|, [l2[ 
and polarization Q as well as techniques for their extraction |l8|, |l9, have been studied extensively in 

the past few years. The four-point function or trispectrum has recently received much attention due to its use in the 
study of the gravitational lensing of the CMB || |2l], |2^] . Techniques for measuring certain components have been 
tested on the Cosmic Background Explorer (COBE) data [^3[ Q . Still, a complete treatment incorporating the full 
symmetry properties of the temperature and polarization fields has been lacking in the literature. 

In this paper we complete the formalism established for the temperature trispectrum plf . The addition of polar- 
ization information leads to a multiplicity of trispectra corresponding to all possible combinations of three observable 
fields. It has been recently shown that the higher point correlations of the CMB polarization contain the majority 
of the information on gravitational lensing in the CMB p5| . Trispectra also quantify the non-Gaussian errors to 
temperature and polarization power spectra measurements. 

The outline of the paper is as follows. We consider the general symmetry and noise properties of trispectra in 
Sec. |l| As an illustration of the construction and noise properties of trispectra, we apply these techniques to the 



initial non-Gaussianity in the curvature fluctuations of the form predicted by slow- roll inflation in Sec. III. We show 
that the sensitivity to initial non-Gaussianity in the trispectra can approach that in the temperature bispectrum [ p"3| . 
In Appendix ^ we summarize relations useful for the study of high order correlations in the polarization. In Appendix 
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H we cover the details in the properties, measurement, and approximation of the trispectra that may be useful for 
future studies. 

II. FORMALISM 

We begin with definitions associated with the harmonic description of the temperature and polarization fields in 



Sec. II A. We conside r the general symmetry properties of n-point harmonic functions in Sec. II B and app ly th em to 



the trispectra in Sec. II C. Finally we derive the Gaussian noise properties of trispectra estimators in Sec. [IE 



A. Definitions 

The temperature field 0(n) = AT(n)/T is decomposed into multipole moments according to 

0(n) = 5>ri7W (1) 

Irn 

The polarization, described by the Stokes parameters Q{n) and U(h) in spherical polar coordinates, is a spin-2 field 
and is similarly decomposed as |^6|, 

(Q±iU)(h) = ^±^r±2^ m (n), (2) 

lm 

where s YJ m (n) is the spin- weighted spherical harmonics whose properties are reviewed in Appendix |a|. Note that 
o*i - H ■ 

Under a parity transformation P taking n to — n, the spin-spherical harmonics transform as 

P[ s Yr(n)} = (-l) l - s Yr(n), (3) 
and so it is convenient to define the parity cigcnfunctions 



ym I y 



* 1 I —s 1 1 

~2i 



s0 rn ^ s I s , _ (5) 



We see that under parity, P[ s £J n ] — (— l) l s £J n , whereas P[ s O™] = (— 1) +1 s O, m . A spin-0 field such as the temperature 
fluctuation carries only the q£ ™ = Y™ mode. The spin-2 polarization field has two components that are distinguished 
by parity M, M 



E?±iBr (6) 



called the E and the B modes. This definition differs from that in |2fJ by an overall minus sign, so in particular, the 
sign of temperature-polarization cross-correlations is reversed. 
The fields E(h) and B(h) on the sky are defined as 



E(h) = J2E7 l Yr(n), (7a) 

lm 

B(n) = Y,BrY l m {n). (7b) 

lm 

The parity properties of the eigenstates s £™ and s O™ imply that E(n) will be a scalar under parity, whereas B(h) 
will be a pseudoscalar. 

Lastly, requiring that the fields 0(n), Q(n), and U(n) be real imposes the constraints 

xf* = (-l) m xr m ,xe {S,E,B} (8) 

on the multipole moments. This constraint also enforces the reality of E(n) and B(fi). 
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B. Rotational and parity invariance 

The harmonic transform of the n-point correlation of CMB fields x, . . . , z € {0, E, B} defines the n-point harmonic 
correlation functions according to 

{x(n 1 )...z{n n )) = ^2{x^...zr n n )Y l ^(h 1 )...Y l ^(h n ). (9) 

Li ,mi 

Since the CMB fields are assumed to be statistically isotropic, the n-point correlations must be invariant under 
rotations. A general rotation R acts on spherical harmonics as 

R\Yr{n)] = Y, Dl mm'(R)Yr'(n). (10) 

m' 

Useful properties of the rotation matrix D l mm , are summarized in Appendix [A|. The n-point harmonic correlation 
must therefore obey 



This invariance demands a specific form for the m-dependence as we shall see. 

Under a parity transformation an n-point function containing k B fields will transform according to 

P[{x™ . . . *£«>] = (-l) k+ Z {x^ . . . (12) 

Invariance under a parity transformation therefore implies that the n-point function containing k B fields will vanish 
when 

k + J2^= odd. (13) 

i 

C. Trispectra 

The reduction of the four-point harmonic function into a rotationally invariant form follows the steps outlined in 



11. Using the Clebsch-Gordan property (A7) on Eq. (fidf) to pair (h,h) and together and applying the 



orthogonality condition (A3) to the resultant pair, we reduce the function to 



WW - £(-D M ( t t 2 —M ) ( i t M ) « <!>). (14) 



LM 



The trispectrum Qy,' 1 ^ (L) represents a configuration with sides l\ . . . I4 labeled by the fields w . . . z, with one diagonal 
of length L forming a triangle with l\ and I2 (Fig. [I]). 

Choosing the two other pairings (Zi, Z3) and (Zi , Z4) yield alternate representations of the trispectra. Since each 
representation is constructed by adding pairs of angular momenta, each representation is complete, and three repre- 
sentations are related through Wigner 6-j symbols via 

Q?S 3 AL) = 5>l)«»-*(2L + l){!j \l L ^Q?S:{L') (15a) 



and 



Q?S: W = E(-i) L+i '( 2i + !) { I I L l } QvS: i L ')- (isb) 

L' 

The trispectrum is obtained by subtracting the unconnected or Gaussian piece from Qy^zi 2 (£), so 



(16) 
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where the Gaussian piece is constructed from the power spectra 



cr = (xrvD 



(17) 



as 



= (-l) il+ V(2ii + l)(2/ 3 + l)C?*C^5 hl2 5 l3ll 8 L0 

+ (2L + 1) [{-l) h+h+L C^Cf:5 hh 5 l2li + C?*C*n hU 5 hh ) ■ 

From the permutation symmetry of the trispectrum (jlj) , additional constraints hold 

= (-1) S -T«(L) = (-l)^TZQ(L) = (-1) E " +E *2^(L), 

where 

S[/ = L + li + I2, 
E L = L + Z 3 + Z 4 . 



(18) 
(19a) 



(19b) 
(19c) 



The constraints ( |l5| ) and ( |l^) express redundancies in the definition of the trispectrum, where a physical configuration 
can be labeled by 4! = 24 different permutations of the field labels and pairings. 

In practice, the following construction guarantees that trispectra obey the constraints outlined above. Given that 
the connected part of the four-point function can be expanded into its three pairings as 



(»rccc)o - EC- 1 

LM 

+ {xi 2 <-> yi 3 ) + {xi 



h h 



L 



mi m,2 —M ) \ rri3 1TI4 M 



h U L 



(20) 



the latter two pairings are projected onto the (h, h) basis through the recoupling relations ( jlq ) to give 



t;S;{l) = p;X !2 (£) + (2l + i)]T 



(-1) 



h h L 
U h L 



l' J h h L 



h U L 



(21) 



Having satisfied the recoupling relations, the remaining constraints (|19|) are enforced by introducing a reduced function 

r m l X! 2 ( L )^ where 



M«3 Z '4 



(l) = ^ ,a ( i ) + (- 1 )^C( i ) + (- 1 )^^C( i ) + (- 1 ) 



(22) 



5 



with the additional constraint 



(23) 



These exhaust the 4! redundancies due to permutations of the fields. Given the functional form of T, the trispectrum 
can be constructed by permuting the fields along with their indices, as indicated above. Therefore T provides the 
most economical description of the trispectra for a given physical effect. 

It is possible to think of the trispectra configurations as being labeled by a fixed field configuration [e.g., Ix^yi' (L) 
with an x and y always related through the diagonal L], with the indices U allowed to vary. The derivation of 
symmetry properties in such a viewpoint is straightforward, and is deferred to Appendix B 1 . The above construction 
automatically enforces the symmetries with fixed field configurations, Eqs. (p3)-(B5). 



D. Noise properties 



By inverting definition (^4|) , we obtain an estimator for a trispectrum 



.M 



h h L 



h U L 



mi 7Tt2 M J \ rvij, f7i4 — M 



K'^fO-^W, (24) 



where the Gaussian estimator G is constructed using expression with the power spectra replaced by their esti- 
mators. We discuss more practical forms of this estimator in Appendix |B 2j . 

The covariance between two trispectrum estimators due to Gaussian noise then becomes 



Sll>n£ + (2L' + I) 



(-1) 



2L+1 

where no two Vs in the primed and unprimed sets are equal, and 



h h L \ AT is,i i\L+l' J h h L 



N- 



:V2 



(25) 



+ [i <-» k,j I], 



(26) 



with fi denoting the field associated with l[. If any two k's are equal, one would need to consider additional terms in 
the covariance arising from pairings within the primed and unprimed sets. 
Using the above covariance, the total signal-to-noise ratio is then given by 



(27) 



li l' L,L' abcdwxyz 



"Cov" here is the covariance in Eq. ( |25| ) viewed as a matrix and the field-type sums are over the measured fields. 

This matrix will possess singular values for permutations that are equivalent. They can be eliminated by singular 
value decomposition or cquivalently by restricting the sums to a set of inequivalent permutations. The latter is 
computationally more efficient and the redundancies expressed in Eqs. ( |l5| ) and ( |f9| ) can be removed by restricting 
the I sums. Thus the total signal-to-noise ratio simplifies to 



v) = E EEEfC^N-'^w), 

h>h>h>h L abed wxyz 
bi 2 



(28) 



where the covariance between T^ 1 ^ 2 (L) and T^ 1 ^ 2 (L) likewise simplifies to 



[Cov] - {2L + \)C%°CgC%C%, 
so that the matrix inverse is now only over field choices for a fixed set of multipoles. 



(29) 
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III. TRISPECTRA FROM INITIAL CONDITIONS 



As an application of the formalism for describing temperature and polarization trispcctra and computing the 
signal-to-noise ratio of their estimators, we consider here the signature of non-Gaussianity that is inherent in the 
initial conditions. In Sees. Ill A and IIIB we motivate a form for the trispectra based on slow-roll inflation. Although 
typical models predict amplitudes far below th at wh ich is potentially observable, this form is generic to local non- 
Gaussianity in the initial conditions. In Sec. Ill C we describe the transfer of this initial non-Gaussianity to the 

HID that the total signal-to-noise ratio in the 



In Sec. 

observable temperature and polarization fields 



We show in Sec 

trispectra is comparable to that in the temperature bispectrum considered previously in the literature 



A. Inflationary motivation 

The standard inflationary paradigm is known to predict a very nearly Gaussian spectrum of initial curvature fluc- 
tuations which under linear gravitational instability theory implies a Gaussian spectrum of the CMB fluctuations. 
However, nonlinear corrections in inflation and gravity can produce non-Gaussian fluctuations, which may be observ- 
able in the microwave background. The imprint of such nonlinearity on the bispectrum of the microwave background 
has been studied extensively. The theoretical predictions for the bispectrum and the related statistic of skewness has 
been described in ||, [l3| [1^], and observational limits placed using existing data in JO], [l8| We here extend 
these considerations to higher order to investigate effects on the trispectra. 

Following , let us consider the non-Gaussianity induced in corrections to the correspondence between a Gaussian 
inflaton fluctuation and the Newtonian curvature $. $ during matter domination can be related to the inflaton 
fluctuation at horizon exit according to 



$0) 



12ttG 



12ttG 



<fio+5<t> 



$0 

dlnff 



d\nH 

d4> 



6ttG d 



dlni? 



2ttG d 2 



dlnH 



d<p 



S4> 3 + 0(84> A ). 



(30) 



The leading order term, which we will denote $^(a;), carries Gaussian random fluctuations from 8<j). The higher order 
terms can then be written as 



(31) 



with 



h 
h 



5 1 d 2 In V 



6 8ttG 

25 1 

54 (8ttG) 2 



d 2 InV 



d 3 In V d In V 

dip 3 d<j) 



(32) 



where V((f>) is the inflaton potential. This model generalizes the considerations of [Q [Tl], [l3|, |2j], 30 1 to higher order. 
Note that our f\ corresponds to f^L in jlq]. As an example consider an inflaton potential of the form V = \<fi n , then 



, 5n fM p 
h = 0, 



(33) 
(34) 



with M 2 = 1/8ttG defining the reduced Planck mass. For inflation to occur, this class of models requires that 
cf> \/120nAfp, so that we obtain /i ~ 0.01 as an order-of-magnitude estimate. 

Although these small coupling coefficients and the observed 10 -5 level of curvature perturbations make non-Gaussian 



contributions from typical inflationary models unmeasurable, the general form in Eq. (31) is simply a perturbative 
expansion of the curvature fluctuations and so we leave /i and ]i as free parameters and explore the extent to which 
they are measurable in temperature and polarization trispectra. Other sources of non-Gaussianity of this form include 
interaction terms in the inflaton potential , stochastic interactions of the long- wavelength inflaton fluctuations with 
the short-wavelength modes || [HJ , and some multiple field models ||| . 
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The ansatz ( |3l| ) for the curvature fluctuations imply higher order correlations in Fourier space since products 
become convolutions. We can decompose the contributions into linear and nonlinear parts, so that 



with 



M*) = A 



where 



l+P 2 +fc), 



(<*>?»> 



P*(fe), 



(35) 

(36) 
(37) 

(38) 



Here Pj>(fc) is the power spectrum of $z,(fc). 



B. Curvature trispectrum 



The connected part of the four-point function (<1? (fci)<t>(A;2 )*E , (fc3)<I>(fe4)} has leading order contributions from terms 
of the form (^A(ki)^A(k 2 )^ L (k 3 )^ L (k 4: )) and ($B(fci)$z,(fc2)$L(fc3)$L(fc4))- To compute the curvature trispec- 
trum, we note that the symmetries with respect to exchange of $(fej) and in the $ four-point function are 
identical to the symmetries with respect to exchange of fields in the CMB trispectrum. We therefore follow the same 
decomposition process by first defining 

($(fc 1 )$(fc 2 )$(fc 3 )$(fc4)) c = (2tt) 3 J d s KS(ki + k 2 + K)S(k 3 + k4-K)T^(ki,k 2 ;k 3 ,k 4 ;K), (39) 

and then constructing 

T 9 (k x ,k 2 ;k 3 ,k 4 ;K) = P # (fc x , fe 2 ; fe 3 , fe 4 ; K) + J d 3 K' [S(k 3 - k 2 - K + K')P 9 (k 1 ,k 3 ;k 2 ,k i ;K') 

+ 5(k 4 - k 2 - K + K')P^(k 1 ,k 4 ;k 3 ,k 2 ;K')}, (40a) 

with P$ constructed out of a reduced trispectrum 7$ according to 

P^{k 1 ,k 2 ;k 3 ,k 4 ;K) = 7i(fci, fc 2 ; fc 3 , fc 4 ; If ) + %>(k 2} k 1 ;k 3 ,k 4 ;K) 

+ T 9 (k 1 ,k 2 ;k i ,k 3 ;K)+T 9 (k2,ki;k A ,k 3 ;K). (40b) 

The leading order contributions to the reduced trispectrum are 

T* A (k u k 2 -k 3 ,k A -K) = 4/?P*(lOP*(fci)P»(*; 3 ), (41a) 
T 4 , r( (fe 1 ,fc 2 ;fc 3 ,fc 4 ;i<r) = / 2 [P$(fc 2 )P $ (fc 3 )P 4 ,(fc4)+P <& (fci)P <& (fc 2 )P <& (fc 4 )]. (41b) 

The curvature trispectrum induces an angular trispectrum onto the CMB fluctuations as we shall now see. 



C. Angular trispectra 



In the linear regime, curvature perturbations generate CMB fluctuations as 

/ii 3 k 
$( k )g al (k)Yr(k), (42) 

where a may be the temperature or E-mode multipolc moment, $(fc) is the primordial curvature perturbation, and 
g a i{k) denotes the radiation transfer function for a = <d,E. The multipolc moments a™ inherit their statistical 
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properties from $(fc), so that in our case, the trispectrum is related directly to integrals of the four-point correlation 
function of <i>(fc). 

From expression ( f42| ) , the harmonic four-point function is related to the $ trispectrum by 



x (2nfg all (k 1 )g b [ 2 (k2)g c i 3 {k 3 )g dl4 (k 4 )T$(k 1 ,k 2 ;k 3 ,k 4 ;K). (43) 

The reduced trispectrum T^ ll d * 2 (L) is then obtained from the reduced $ trispectrum simply by replacing T$ in the 
above relation and performing the integrals over directions ki and K, so that 

Kit! W = (~) f rldnrldn (kldh) ■ ■ ■ (kldki) K 2 dKj L (K ri )j L (Kr 2 ) 

X [gah( k l)jh( k l r l)} [9bl 2 ( k 2)jl 2 (k2ri)} [9cl 3 ( k 3)jh( k 3 r 2)} [ff<ii 4 ( k i)jl 4 ( k 4r 2 )} 

x T^(ki,k2;k3,k4,;K)hi 1 u 2 hi 3 LU, (44) 



where 



/ (2^ 1 + l)(2l 2 + l)(2L + l) ( h l 2 L 
MlLl2 ~ V 4vr ^000 

and ji{x) are the spherical Bessel functions. 

Substituting expressions ( pl| ) into the above, we find that the reduced trispectrum is given by 



(45) 



with 



Kit (L) = T A %fc (L) + T B %% (L), (46a) 



T ^S L ) = y^drx^^Fi^.raJaf^ri^^Jof,^)/^^)^^,,^, (46b) 
TbZI: (L) = f r 2 d<(r)/3f 4 (r) (r) + /9£ (r)^ 3 W] Kl^u, (46c) 



and 



i^(ri,r 2 ) = 




/" K 2 dKP*{K)j L (K ri )j L (Kr 2 


«f(r) = 




[ fc 2 dfc(2/ 1 ) 5ai (fc) Ji (fcr), 


/TO = 




/" kHkP^{k)g al {k)ji{kr), 






f k 2 dkf 2 gai( k )ji(kr). 



(47a) 
(47b) 
(47c) 
(47d) 

The trispectrum is formed using Eqs. @ and (§§). 

To properly evaluate the trispectra, one must extract the radiation transfer function gi{k) numerically from an 
Einstein-Boltzmann solver as has been done for the temperature bispectrum Jl3| . This process is numerically cum- 
bersome and we instead seek an analytic approximation of its effects. 

For small multipole moments (I <C 100), CMB temperature fluctuations arise mainly from the Sachs- Wolfe effect 
||| , Here the radiation transfer function g@i(k) takes on the simple form 

gei(k) = pi(kr.), (48) 

with r* = r/o —ri xec denoting the conformal time elapsed between recombination and the present. In this regime, otf (r) 
and nf(r) simplify to 

a ?<r) = ^ir-r,), ( 49 ) 



m?W = ( 5 °) 
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FIG. 2: Signal-to-noise ratio in the temperature trispectrum as a function of the maximum multipole Z max in the Sachs- Wolfe 
approximation (Eq. ^). The dotted line corresponds to /i = = 0, the dashed line to fi = Q,f 2 = 1, and the solid line to 

/l = l,/2 = l. 



Since the temperature power spectrum is given by 

Cf W = k 2 dkPt,(k)j?(kr*), (51) 

the other functions can be related to Cf w as -Fl (r* , r* ) = 9Cf w and (3f(r*) = 3Cf w . The reduced trispectrum can 
then be expressed in terms of Cf w as 

(L) = 9Q s 2 w Cf 4 w [ 4/l 2 d w + h K w + C™)] h h u 2 h hLh . (52) 

Signal-to- noise ratios calculated from expression (|5^) for various choices of /i and are shown in Fig. (||). 

To estimate the trispectrum for higher multipoles note that for the /i term in (|52|), Cf 2 w and C^ w appear from an 
integration over radiation transfer functions that is very similar in form to that of the true power spectrum C ; ee . An 
approximation to the trispectrum induced by /i then becomes 

T^(L) » 36h hLl2 h hLl j!C^C^C^ (53) 

and should be valid to the extent to which the anisotropies result from slowly varying local temperature fluctuations on 
a thin last scattering surface. While the /2 term does not have this simple form, we take the /i piece as representative. 

The Newtonian curvature $(fc) also acts as a source for E'-mode polarization, through the anisotropy of Compton 
scattering which links the local quadrupoles of temperature fluctuations to local -E-mode fluctuations. Although there 
is no equivalent to the Sachs- Wolfe approximation for the E-mode radiation transfer function, we can again take the 
above approximation for a slowly-varying source to obtain 

T A E E f E \l(L) » 36h tlUa h hLU ftdF r CZ E (% B i (54) 

where C[ w is still the Sachs- Wolfe approximation to the temperature power spectrum. Thus, for /2 = 0, the _E-mode 
trispectrum should behave similarly to the temperature trispectrum. Mixed and E trispectra would take on an 
analogous form. 



D. Signal-to-Noise 



We utilize the formalism described in Sec. IIP to calculate the expected signal-to-noise ratio for primordial non- 
Gaussianity. Setting /2 = 0, we use Eq. ([33) for the temperature trispectrum and compute (S/N) 2 to the cosmic 
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10 100 1000 
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FIG. 3: (S/N) 2 / ft vs. /max in the temperature trispectrum, using the approximation (|53|). 

variance limit, with Gaussian contamination from gravitational lensing. Figure ^ shows (S/N) 2 / ff as a function of 
the maximum multipole moment observed l max . The tapering of (S/N) 2 is due to the fact that the noise contribution 
from lensing becomes significant at small angular scales. The figure indicates that the temperature trispectrum may 
be sensitive to non-Gaussianity for fx < 10, although a detailed calculation using expressions (|4^) involving the full 
radiation transfer function will be necessary to place rigorous bounds. Since polarization trispectra take a form similar 
to the temperature, (S/N) 2 can be at most enhanced by the number of independent trispectra terms. For E and 
combinations this represents a factor of a few. 

Compared with the sensitivity of the temperature bispectrum to primordial non-Gaussianity, the temperature and 
polarization trispectra contain a comparable amount of information The (S/N) 2 in the trispectra can exceed 

that of the bispectrum if fx 1 due to the steep scaling of ff. 

IV. DISCUSSION 

We have introduced a complete formalism for the study of 4-point correlations in the CMB temperature and 
polarization fields. This formalism should be useful in future tests of the non-Gaussianity of the CMB induced in the 
early universe and by the evolution of structure. It is also of use in determining the non-Gaussian contributions to 
errors in temperature-polarization power spectra measurements. 

We have applied these techniques to a particular form of trispectra motivated by inflation, generalizing previous 
treatments to higher order in the initial nonlinearity of the curvature fluctuations. Typical slow-roll inflationary models 
predict an amplitude to the trispectra that is far from observable and so a detection of this type of non-Gaussianity 
would rule out a large class of models. We have shown that because of the large number of trispectra configurations, 
the sensitivity to initial non-Gaussianity in the trispectra approach that of the well-studied temperature bispectrum 
at high multipoles. 

Trispectra from secondary anisotropies such as gravitational lensing |2lJ are expected to be substantially larger 
and should be fruitful ground for future studies. While measurement of these non-Gaussian signatures will no doubt 
prove challenging due to foregrounds, systematic effects and computational cost, the wealth of information potentially 
contained therein may justify the large effort that will be required. 
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APPENDIX A: SPIN- WEIGHTED SPHERICAL HARMONICS 

We summarize the conventions and properties related to the spin- weighted spherical harmonics p3]. A function 
s f(0<f)) on the sphere carries a spin weight s if, under a right-handed rotation of the basis (eg, e^) by an angle ip, it 
transforms according to s f(0(j)) — > e~ %s ^ ' s f(0<p). For such functions, there exist a complete and orthonormal basis with 
spin weight s, called the spin- weighted spherical harmonics. These spin-s spherical harmonics s Y[ m can be constructed 
from the ordinary spherical harmonics by application of the raising and lowering operators |34f| . Alternately, they are 
given in terms of rotation matrices as 



s Yr(M = (-iy^^±lD l _ ms (a(3 7 )e^. (Al) 

The Euler angles specify a rotation around the coordinate z axis by 7, followed by a rotation by (3 about y, then a 
rotation by a about the (original) z axis. The rotation matrix is given explicitly by (see, e.g., |3q]) 



Dl m8 (<*h) = 



(l + m)\(l-m)\ 



sin 2i (/3/2) 



x 

k 



(i + s y.(i- s y. 

T( l V)( k i:: m )(-^ +l+s ^ + -(m- (A2) 



Note that this convention for S Y" 1 differs from that presented in J34|,j37j by (— l) m , but corresponds to the Condon- 
Shortley convention for the ordinary spherical harmonics when s = p6[| . Below and throughout the paper we use the 
shorthand convention for the arguments of the spin-spherical harmonics and rotation matrices ft = (6<f>), R = (a/37) 
and their differential elements dft = defy d cos 9, dR = da d cos /3dj. 

The properties of the spin-weighted spherical harmonics follow from those of the rotation matrices. In the text, we 
utilize four such properties: orthogonality, completeness, Clebsch-Gordan expansion, and angle addition. 
Orthogonality: 



implies 



Completeness: 



implies 



J dRD l ; ns (R)D l m , s ,(R) = -^jS w S mm ,5 ss , (A3) 



dn s ^*(n) s ^ (n) = S lv S mm ,. (A4) 



^D£ s (a/3 7 )£>L(aW) = ^^6(a - a')6(co S p - cos 0')5(~f - j') (A5) 

hns 



J2sY l m *(9(f,) s Y l m (9'(t)') = 5(<f>-(j>')6(cos9-cos9'). (A6) 



lm 

Clebs eh- Gordon relation: 

h h L \ ( h h L \ , _, ,m+m' t^l 



LMM' 

implies 



D h ,(R)D l > ,(R) = Y (2L+1) 11 12 \. )[ ll , l2 , \„ )(-l) M+M 'Db M ,(R), (A7) 



sl Y^(n) S2 Y lT (n) = VWi + 1)(2Z 2 + 1) £ (~l) M+S { _^ _^ £ J ( £ J _ 5 J ^(n) <A8) 

Auxiliary (orthogonality-Clebsch-Gordon) relation: 

[dRD h ,(R)D h ,(R)D h ,(R) = 8tt 2 ( h h h )( h , h , h ,), (A9) 



implies 



dn sl Y^(n) S2 Y™(n) S3 Y™°(n) 



4tt 



l[ 2k + 1 



1/2 



'i h h 



S 1 -S 2 -S3 / V TO 1 TO 2 TO 3 



^1 ^2 ^3 



if sx + S2 + S3 = 0. Addition theorem: 



implies 



2/ + 1 

47T 



(-l)^ S2 ^- Sl (/3a)e is ^ 
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(A10) 



(All) 



(A12) 



The relationship between the angles (shown in figure is given explicitly by 

sin0' 



cot a 



cos 8' cotU' - <f) + cot 9— 



cos (3 = cos cos 0' + sin0sin#'cos((/>' 



cot 7 = cos 9 cot (4>' — 4>) — cot 6 



sin ( 



sin 



(A13) 



The addition relation corrects a sign error in |37fl and agrees with p8j, once one accounts for the differences in the 




FIG. 4: Relation between Euler angles (a/37) an< i the original rotation angles (9(f)) and (6' '<j>') for the weighted sky maps (B8), 
with the identification n — > (9' <f>'), q — ► (00), so that 0^ = a, and n ■ q = f3. 

phase convention for glT™. 



APPENDIX B: ADDITIONAL TRISPECTRA PROPERTIES 



1. Symmetries 



Here, we present the symmetry properties of trispectra when only the angular momentum labels are permuted, 
keeping the field labels fixed. Such a representation is redundant but can be useful if the diagonal length L is related 
to some physical quantity, as is true for the inflationary trispectra [see Eq. (|4^)]. 
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Using the three CMB fields {8, E, B}, 15 distinct four-point functions can be constructed, with the field contents 

xxxx e {0969, EEEE, BBBB}, 

xxxy e {QQQE, 999B, EEEQ, EEEB, BBBQ, BBBE}, 
xxyy e {QQEE, QQBB, EEBB}, and 

xxyz e {QQEB, EEQB, BBQE}. (Bl) 

For each case, the use of identical fields results in a restricted set of permutation symmetries. 

The case for (9990) has been worked out in detail in pl| , where the trispectrum was shown to be composed 
from a reduced form, thus incorporating the 4! = 24 possible permutations of Zj's which leave the four-point harmonic 
function unchanged. Following the treatment in pi] , f° r trispectra of the form xxxx, permutation symmetry of ij 
requires that the trispectrum Qx',lx',l(L) = obey the constraints 

Q\±{L) = (-l) s ^Q^(L) = (-l) Si Q^(i) = (-l)^+^Q^(£), (B2a) 

Ql±( L ) = E(- 1 ) /3+/3 ( 2i + 1 ){!4 \\ i}Q l &S L ')> ( B2b ) 



and 

= ^(-l) L+L '(2i + l)|;; g ^ ^Q^(i'): (B2c) 

where £;y = i + l± + I2, and £l = L + Z3 + £4. 

For trispectra of the form xxxy, the permutations (123) on the Z indices are allowed, so that the constraints are 
given by 

Ql'SliL) = {-l)*"Ql\lll{L) (B3a) 
= X>l) i2+(3 (2L + 1) j \\ \\ V }^v!:(i') (B3b) 



^(-l) L+i (2L + 1) I g g J I Q*!;^ (£'), (B3c) 

where the last two relations come from Eq. (|l^). These exhaust the allowed 3! = 6 permutation symmetries. 

For trispectra related to (x™ 1 y™ 2 x™ 3 y™ 4 } , the permutations (7i <-> Z 3 ) and (Z2 *-* U) are separately allowed, so that 
there should be 2 • 2 = 4 permutations to account for. In the xxyy basis, the permutation symmetries imply that 

Qyllml (L) = (-lf u Q7 3 yll (L) = (-1) E 'C ( L ) = ("l^C ( B4 ) 

Expressions in the other bases can be derived through the use of the recoupling relations ( |l5| ) . 

Lastly, for trispectra related to four-point functions (x™ 1 x™ 2 V™ 3 z™ 4 ) the only permutation symmetry allowed is 
the exchange of l\ and 1%. Accordingly, the trispectra obey 

QZ&W = (-1)*°Q%%(L), (B5) 

where, again, the recoupling relations can be used to find the symmetry constraints in the other bases. 

The signal-to-noise ratio for trispectra with fixed field configurations can be obtained from expression ( pfil ) by using 
the recoupling relations ( fl5| ) to permute the field symbols into that of the fixed field representation. The results can 
be simplified by using the following identities for the 6-j symbols, 

^(2e + i)l:5 e f U:5:)-^ m 



4" I c d f J I c d 9 / 2/+1' 

and 



E(-')- +/+ ^+'){:s;}{s^}={^^ 
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Rewriting the signal-to-noise ratio in terms of the fixed field configurations transfers the redundancies in the field 
permutations into redundancies in the I configurations, so that the restrictions on the sum over l^s typically become 
relaxed. 

2. Measurement 

Direct measurement of the trispectrum using the estimator of Eq. is computationally expensive due in part to 
the quintuple sum over m's. Since the m dependence of the four-point function simply reflects the rotational invariance, 
it is useful to find an estimator that employs these symmetries in a more efficient way. The following construction 
parallels that of the temperature bispectrum jl7) and trispectrum ||l| and takes into account the subtleties due to 
the spin-2 behavior of the polarization fields. For the temperature and polarization fields, one can define a set of 
weighted sky maps ef(q), where 



e?(q) = / ,hiU[n)IMn-q). (BSa) 



and 



ef(<7) = y 4^ [1 + 2)1 / dflP ii h ' <?) [Q(n) cos 2<^ + U(h) sin 20 ft ] , (B8b) 
1 21 + 1 (I — 21! C 

ef (q) = J 4?r j + 2 jj J dhP?(h ■ q) [U(h) cos 2^ - Q(n) sin 2fo] . (B8c) 

The angle tpa is the angle between the great circles defined by (n, z) and (n, q), and serves to transform the Stokes 
parameters from the spherical polar basis to the great circle basis (see Fig. [|). Using these maps, the quantity 

q7Sii (£) = tZS: w + g7S: m m 

can be estimated by expanding the harm onic coefficients in the direct estimator ( ^4|) back into fields, expanding the 



Wigner 3-j symbols using equation (A10), resulting in the expression 

(S o) (o q)G%%( L ) = ( 2L + 1 ) / ^^^(n 1 -n 2 )e^(n 1 )ef 2 (n 1 ) e f 3 (n 2 )ef 4 (n 2 ). (BIO) 
The 3-j symbols impose the constraint ii + i 2 + L =even, so the above expression does not allow measurem ent o f 



modes with l\ + 1% + L =odd. It may be computationally advantageous to compute the double integral in Eq. (B1C) 
as a single sum in multipole space by noting that they individually return the harmonic decomposition of a product of 
two e-fields p3f . One must also account for complications arising from realistic issues, including the leakage between 
the E and B modes due to incomplete sky coverage, for example, by Monte-Carlo techniques. 



3. Flat sky approximation 

A sufficiently small patch of sky can be considered flat. In this limit it is computationally and conceptually 

advantageous to consider the Fourier representation of the trispectrum. Here wc establish the relationship between 
the angular and flat-sky trispectra. 

In the flat-sky approximation, the temperature and polarization fields are expanded in Fourier modes as 

e(n) = J -^Le(l)e il -*, (Blla) 

±A(h) = - J ^L ± A(l)e ±2 ^^e a - A , (Bllb) 

where 4>i is the azimuthal angle of I, and the Stokes parameters are defined in a spherical basis. Again, the E and B 
modes are defined as 

±A(l) = E(l)±iB(l), (B12) 
and the two-point angular correlation functions reduce to 

(x*(l)x'(l>)) = (2nfS(l-nCf l f. (B13) 
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The connected part of the four-point correlation function can be written as 

(w(h)x(l 2 )y(l 3 )z(h)) c = (2ir) 2 S(l 1+ l 2 +l 3 +l A )T^ l 3 Q\l 12 ,l 13 ), (B14) 

where h 2 and h 3 denote the lengths of the two diagonals. This can be broken up into pieces corresponding to distinct 
pairings, so that 

Tt:;: A )\hM = <^fa 12 )+<3 ) a 1 3)+p£;s ) ^)> ( B15 ) 

with Z14 being a function of h 2 and h 3 . The second and third terms can be projected onto the first pairing. Denoting 
the trispectrum with the projected terms as T/™ 11 *' 2 ? (L), the four-point function becomes 

(w(h)x(l 2 )y(l 3 )z(h)) c = (2ir) 2 J d 2 LS(h + l 2 + L)6(l 3 +l 4 - L)T { ^^{L), (B16) 
where we used a decomposition of the delta function 

S(h +l 2 +l 3 + h) = J d 2 LS(h +l 2 + L)S(l 3 +h- L). (B17) 



Separating the pairings according to the prescription in Eq. (B15), the flat-sky four-point function is separated into 



(w(h)x(l 2 )y(l 3 )z(l 4 )) c = (2tt) 2 J d 2 L {d(h +l 2 + L)6(l 3 + h~ L) P^f (L) 
+ 5(h + l 3 + L)S(h +h- L)P^^(L) 

+ 8{h + h + L)5(h+h-L)P { ™^ (L)}. (B18) 

To relate the above expression for the four-point function in the flat sky approximation to the full-sky expression, 
we use the relation between the Fourier coefficients x(l) and x™, 



= d '^i;^^ (B19a) 

m 

and the inverse 

-r m / ^- e - tm<l "x(l), (B19b) 

J 27T 



21 + 1 

An 



where x £ {0, E, B}. The phase convention chosen here differs from those in ^6), due to differences in the choice 
of phase in the definition of the spin-weighted spherical harmonics. 

The plane waves in the 5 functions are further decomposed into spherical harmonics using the relation 

2 j J2{-i) m Yre- m<t " , (B20) 

m 

where the approximation is valid for small angles (or large multipoles I). 

If the four-point function has an even net parity, then Py^zi 2 (L) is independent of the orientation of the quadri- 
laterals, so that the integrals over the azimuthal angles 4>i can be performed, and we obtain the desired relation 

« 2 W « ^/(21 1 + D...(2Z 4 + 1) ( £ l £ L ) ( ^ L Q ) PfctfiL). (B21) 

For odd net parity, the Wigner 3-j symbols should be reinterpreted as their analytic continuation (see Q), Eq. (Bl) 
and Appendix C3). 
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